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ABSTRACT 

We present B-, V-, R^-, Ic-, J-, H-, K- and K -band observations of the optical transient (OT) associated with 
GRB 970228, based on a reanalysis of previously used images and unpublished data. In order to minimize cal- 
ibration differences we have collected and analyzed most of the photometry and consistently determined the 
magnitude of the OT relative to a set of secondary field stars. We confirm our earlier finding that the early decay 
of the light curves (before March 6, 1997) was faster than that at intermediate times (between March 6 and April 
7, 1997). At late times the light curves resume a fast decay (after April 7, 1997). The early-time observations of 
GRB 970228 are consistent with relativistic blast-wave models but the intermediate- and late-time observations 
are hard to understand in this framework. The observations are well explained by an initial power law decay 
with a = -1. 73^0 [2 modified at later times by a type-Ic supernova light curve. Together with the evidence for 
GRB 980326 and GRB 980425 this gives further support for the idea that at least some GRBs are associated with 
a possibly rare type of supernova. 

Subject headings: gamma rays: bursts - supernovae: general 



L INTRODUCTION 

The 7-ray burst of February 28, 1997, was the first for which 
a fading X-ray (Costa et al. 1997) and an optical counterpart 
(Groot et al. 1997a; van Paradijs et al. 1997) were found. 
After the counterpart had weakened by several magnitudes, it 
was found to coincide with an extended object (van Paradijs et 
al. 1997; Groot et al. 1997b; Metzger et al. 1997a). In sub- 
sequent observations with the Hubble Space Telescope (HST) 
on March 26 and April 7, 1997, it was found that the opti- 
cal counterpart consisted of a point source and an extended 
(~ 1") object (GAL), offset from the point source by ~ 0.5" 
(Sahu et al. 1997). The extended object is most likely the host 
galaxy of GRB 970228 and its redshift was recently determined 
at z = 0.695 ± 0.002 (Djorgovski et al. 1999). The original 
light curve was compiled from data taken from the literature 
(Galama et al. 1997). Data obtained until April 1997 suggested 
a slowing down of the decay of the optical brightness. However, 
the HST observations of September 1997 (Fruchter et al. 1997, 
1999; Castander and Lamb 1998a) showed that after April 1997 
the light curve of the point source roughly continued a power- 
law decay. 

We have collected and reanalyzed the photometric informa- 
tion (previously drawn from the literature and presented by 



Galama et al. 1997), collected yet unpublished images and 
present here reassessed Cousins V-, Rc-, and I^-band optical 
light curves and spectral flux distributions of the GRB counter- 
part (Sect, p] , ^,]3]and Q). We confirm the result of Galama et 
al. (1997) that the early decay of the light curves (before March 
6, 1997) was faster than that at later times (between March 6 
and April 7, 1997) and discuss the observations in terms of the 
relativistic blast-wave and blast-wave plus underlying super- 
nova models in Sect. ^. We find that the observations are well 
explained by an initial power law decay with a = -l-73;!;o j| 
modified at later times by a type-Ic supernova light curve. This 
is in agreement with a recent suggestion by Reichart (1999). 
Together with similar evidence for GRB 980326 (Bloom et al. 
1999) and GRB 980425 (Galama et al. 1998) this gives further 
support for the idea that at least some GRBs are associated with 
a possibly rare type of supernova. 

2. REANALYSIS OF PREVIOUS OBSERVATIONS 

We collected and analyzed most of the photometry of 
GRB 970228. In Table |l| we provide a log of all the photom- 
etry on GRB 970228, obtained in the B, V, Rc, and Ic pass- 
bands (corresponding to the Johnson B and Cousins VRI sys- 
tem), some photometry derived from unfiltered passbands and 
translated to B or Rc, the HST WFPC2 and STIS filters, and 



Astronomical Institute 'Anton Pannekoek', University of Amsterdam, & Center for Higli Energy Astrophysics, Kruislaan 403, 1098 SJ Amsterdam, The Nether- 
lands 

^Department of Physical Sciences, University of Hertfordshire, College Lane, Hatfield, Herts, ALIO 9AB, UK 
^Department of Physics and Astronomy, SUNY Stony Brook, NY 1 1794-3800, USA 
^Physics Department, University of Alabama in Huntsville, Huntsville AL 35899, USA 
^Universities Space Research Asociation 
''NASA/MSEC, Code SD-40, Huntsville AL 35812, USA 

'Space Telescope Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA 
*CNR Bologna, Via R Gobetti 101, 40129 Bologna, Italy 

'Copenhagen University Observatory, Juliane Maries Vej 30, D-2100, Copenhagen A, Denmark 
'"Astronomy Dept., University of Washington, Box 351580, Seattle WA 98195-1580, USA 
"SIRTE Science Center, Caltech, MS 314-6 Pasadena CA 91 125 

'-Thiiringer Landesstemwarte Tautenburg, Stemwarte 5, D-07778, Tautenburg, Germany 

1 



2 



the near-infrared J, H, K and K passbands. In order to min- 
imize calibration differences as a result of use of different fil- 
ters/instruments and differences in analysis we consistently de- 
termined the magnitude of the optical transient (OT) relative to 
our calibration, which is described in Sect. 2.1, using a set of 



secondary field stars (Table ||). In Sect. 2.2 and 2.3 we discuss 
each observation separately. 

2.1. Calibration 

As part of the ESO Very Large Telescope Unit 1 (VLT- 
UTl) Science Verification (1998) a series of deep images of the 
GRB 970228 field were obtained in B, V and Rc (see Table []). 
The images were bias-subtracted and flat-fielded in the standard 
fashion. Through the entire night observations of the Landolt 
fields SA 99 (around star 253), PG1633H-099, MARK A and 
TPHE (Landolt 1992) were obtained allowing accurate photo- 
metric calibration for a number of reference stars (see Table ^. 
We determined the transformation of the reference stars' mag- 
nitude to the Landolt system by fitting for the zero point and the 
first order extinction and color coefficients. 

We also determined a calibration in B, V and Rc from 
observations of the Landolt fields SA98 (around star 978), 
PG0918H-029 and Rubin 149, taken as part of the Isaac Newton 
Telescope (INT) observations of March 8 and March 9, 1997 
(see Table |l]). These two calibrations agree, for stars C-H, to 
within 0.04, 0.02 and 0.10 magnitudes for B, V and Rc, respec- 
tively. The difference between the two calibrations is due to 
lesser depth of the INT calibration and we estimated that the ab- 
solute calibration uncertainty in each filter from the VLT-UTl 
observations is 0.05 magnitudes. We used the VLT-UTl mag- 
nitudes for the reference stars C-H and included stars A and B 
from the INT calibration (stars A and B are not contained in 
the field of the VLT-UTl images), correcting for the small zero 
point offset between the VLT-UTl and INT calibration. The 
adopted reference stars' magnitudes are listed in Table ^ 

Observations of the field of GRB 970228 in B, V, Rc and Ic 
(for 240 sec each filter) and of Landolt field SA 98 were ob- 
tained March 2, 1999 with the 1-m Johannes Kapteyn Telescope 
(JKT; La Palma). Of stars A-H only star B was well detected, 
because of the lesser depth of the JKT images. Therefore, to 
obtain the Ic-band calibration, we proceeded as follows. First, 
we determined the magnitudes of star B and six more bright 
stars in the field of GRB 970228, correcting for the zero point, 
the atmospheric extinction and also for a first-order color term. 
We then used the Ic-band calibration of these bright field stars 
to calibrate stars A-H using the William Herschel Telescope 
(WHT) March 8, 1997 Ic-band image; these magnitudes are 
listed in Table ^ For consistency: the calibration in B, V, and 
Rc for star B is in good agreement with the calibration deter- 
mined from the VLT-UTl and INT data (better than 0.07 mag). 

Near-infrared images in the J and K band of the field of 
GRB 970228 and of standard star FS 14 (Casali and Hawarden 
1992) were obtained on March 30.2, 1997 UT with the Near 
Infrared Camera (NIRC) at the Keck-I 10-m telescope. From 
these observations we obtained a calibration of the nearby star 
(starH); we found 7 = 20.27 ±0.10 and 19.67 ±0.10. An- 
other set of near-infrared images of the field of GRB 970228 
and of standard star HD 84800 (Casali and Hawarden 1992) 
were obtained in J, H and K with the Calar Alto 3.5-m tele- 
scope on March 17.8, 1997 UT. From these observations we 
found for the nearby star (star H): J = 20.20 ±0. 15, // = 19. 16 ± 

'^Note that in Galama et al. 1997 the time of the observation of Guamieri et al. 



0.15 and K = 19.59 ±0.28. The two calibrations are in good 
agreement with each other and we adopted their weighted aver- 
ages: 7 = 20.25±0.10,//= 19. 16 ±0.15 and 19.66±0.10. 

2.2. The details of each observation 

The earliest image of the OT was reported by Pedichini et 
al. (1997). The reported magnitude was given relative to the 
nearby late-type star (LTS; star H in Table 0). This observation 
was obtained without a filter; we have transformed this unfil- 
tered magnitude to the Rc band, using the reported filter char- 
acteristics (Pedichini et al. 1997), taking V-Zc = 0.50 ± 0.23 
(from the February 28.99, 1997 UT WHT observations) and 
assuming a power-law spectrum ccv^ . The resulting value 
/?c = 20.0 ± 0.5 is consistent with that reported by Pedichini et 
al. (1997) and Galama et al. (1997). 

In the images of Guarnieri et al. (1997)'-', the OT was 
blended with the nearby late-type star We analyzed their 
February 28.83, 1997 UT Rc-band image and used a large aper- 
ture to include both the OT and the late-type star. The resulting 
magnitude = 20.30 ± 0.28 (OT + GAL + LTS) in our calibra- 
tion is about 0.8 mag brighter than that reported by Guarnieri 
et al. (1997); correcting for the LTS we found = 20.58 ± 
0.28 (OT + GAL). We have not included the subsequent mea- 
surements and upper-limits as given by Guarnieri et al. (1997) 
as they are consistent with detections of the late-type star only. 

Reanalysis of our WHT, INT and New Technology Telescope 
(NTT; La Silla) images have led to marginal corrections; all 
corrections, apart from the V band February 28, 1997 (which is 
0.3 mag off), are within the errors as presented in Galama et al. 
(1997). We also present an unpublished Rc-band non-detection 
(INT March 8.89, 1997 UT). 

TheMarch3.1, 1997 UT observation by Margonetal. (1997) 
with the Astrophysical Research Consortium's 3.5-m telescope 
at the Apache Point Observatory was performed using an un- 
filtered (UF) passband (blue). As the OT was faint in the APO 
image we determined its position relative to a number of nearby 
field stars in the February 28.99, 1997 UT V-band WHT image 
(where the OT is bright) and performed a rotation plus trans- 
lation to obtain the corresponding pixel position in the APO 
image; subsequently we have used this position to center the 
aperture to determine the OT's magnitude (the aperture was 
chosen such that the extended emission was also contained). 
The unfiltered OT's magnitude was translated to the B pass- 
band by fitting an empirical linear relation between the UF -B 
magnitude difference and the B-V magnitude difference de- 
termined from the reference stars A-H (Table We found 
UF-B = 24.03(65) + 0.45(1 l)(B-y) (x^ = 1.2 with 3 d.o.f). 
We obtained a rough estimate of the color of the OT by fitting a 
temporal power law decay, /v oc f"", to the V-band data, inter- 
polating the fit to March 9.85, 1997 UT and using the B-band 
measurement at that same epoch; we found B-V = 1.10±0.34 
on March 9.85, 1997 UT. Assuming that the color of the OT 
did not substantially change this corresponds to Z? = 24.5 ± 0.7 
on March 3.1, 1997 UT. (Note that the error is not dominated 
by the assumed color). This is about 1.4 different from the 
original value (Margon et al. 1997). 

Apart from the March 4.9, 1997 UT V-band image (presented 
by Galama et al. 1997) additional, and yet unpublished, unfil- 
tered (red) and B-band images were obtained with the Nordic 
Optical Telescope (NOT). The OT was clearly detected in the 
unfiltered March 4.9, 1997 UT image but not in the V and B- 

1997 is given incorrectly. 
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band images of March 4.9, 1997 UT nor in the unfiltered im- 
age of March 3.9, 1997 UT. Similar to the procedure for the 
APO data we determined the pixel position of the OT from the 
February 28.99 V-band WHT image to obtain the magnitude of 
the OT in the unfiltered NOT March 4.9, 1997 UT image. The 
unfiltered OT's magnitude was translated to the Rc passband by 
fitting an empirical linear relation between UF-Rc and V -Rc 
determined from the reference stars A-1 (Table ||). We found R^ 
= 22.97(26) + 0.4Eil6)iV-Rc) (x^ = 4.3 with 6 d.o.f). Taking 
V-Rc = 1.03 (corresponding toV-Ic = 2.24 ± 0.14, as ob- 
served for the OT with the HST March 26 and April 7, 1997; 
see below) we found R^ = 23.46 ± 0.32. In a similar way, for 
V-Ic = 0.50 ± 0.23 (WHT February 28.99, 1997) we obtained 
Rc = 23.10 ± 0.27; we adopted the mean value Rc = 23.28 ± 
0.40. Using the same empirical relation between UF-Rc and 
V-Rc we obtained > 22.7±0.2 for the unfiltered March 3.9, 
1997 UT image (using a conservative V - /c = 2.24). 

We also determined the pixel position of the OT to obtain its 
magnitude in the Keck images of March 6.3 and April 5.8, 1997 
UT, using again the February 28.99, 1997 UT V-band WHT im- 
age. We found the March 6.3 Rc magnitude, relative to our cal- 
ibration, to be 0.4 mag brighter than that reported by Metzger 
et al. (1997a); the April 5.8 magnitude is consistent with that 
reported by Metzger et al. (1997b). 

We have not analyzed the Palomar 200 September 4, 1997 
observation of Djorgovski et al. (1997); it is consistent with 
detection of the nebulosity only. 

We reanalyzed the GRB 970228 HST WFPC2 (F606W and 
F814W) and STIS images of Fruchter et al. (1999). For the 
details of the observations we refer the reader to that work. We 
determined the WFPC2 magnitudes of the OT and stars C, G 
and H, using a 2 pixel aperture radius. We corrected the in- 
strumental magnitudes for: (i) the 'long versus short anomaly' 
(Casertano and Mutchler 1998; corrections are -0.16, -0.24, - 
0.12 and -0.19 mag for the WFPC2 observations of March 26 
(V), April 7 (V), March 26 (1), April 7 (1), respectively), (ii) the 
effect of non-optimal Charge Transfer Efficiency (CTE; Whit- 
more and Heyer, 1997; Whitmore 1998; corrections are ~ - 
0.06 mag), and (iii) the geometric distortion of WFPC2 images 
(Holtzman et al. 1995a; corrections are smaller than 0.03 mag). 
Next we determined the aperture correction from 2 to 1 1 pix- 
els, used the instrumental zero points from Baggett et al. (1997) 
and the color corrections from Holtzman et al. (1995b) to ob- 
tain the final V- and Ic-band magnitudes. The errors were deter- 
mined by adding in quadrature the measurement error (Poisson 
noise), the errors in the above-mentioned corrections, the aper- 
ture correction error and a 2% uncertainty for the instrumental 
zero points. For the OT the error in the magnitude is dominated 
by the 'long versus short anomaly' correction (the correction 
is larger for fainter objects). The results for the OT are given 
in Table |l|. We found systematically brighter magnitudes than 
Fruchter et al. (1999) and Castander and Lamb (1998a), which 
can be explained by the 'long versus short anomaly' correction 
that we applied. We found the OT to be redder (V-Ic = 2.22 ± 
0.25 and y-/c = 2.26 ± 0.30, March 26 and April 7, 1997 UT, 
respectively) than found by Fruchter et al. (1999), but consis- 
tent with the result of Castander and Lamb (1998a). For stars 
C, G and H we obtained (V, Ic) = (21.76 ± 0.05, 20.10 ± 0.05), 
(23.46 ± 0.05, 22.52 ± 0.07) and (22.85 ± 0.04, 21.04 ± 0.04), 
respectively. Comparing this with the ground-based calibration 



(see Table g) we found the agreement for the bright star C to be 
excellent (better than 0.02 mag), for star H we found 0.15 and 
0. 12 mag differences and the V-band HST magnitude of star G 
was 0.22 mag off. 

The HST STIS observations were taken in unfiltered mode, 
meaning that the response covers a wide wavelength range from 
about 2000 to 10000 A. To convert the observed counts into 
broadband magnitudes we had to assume therefore something 
about the spectral energy distribution. Our procedure was as 
follows: (i) measure the counts in a 2-pixel radius aperture; 
(ii) add an aperture correction of -1.15 mags based on brighter 
stars in the field, assuming that 10% of the light falls outside a 
large 1 1 pixel radius aperture; (iii) assume a power law slope 
with spectral index chosen to be consistent with the earlier 
WFPC2 observations (/3 = -4.7), and use the SYNPHOT rou- 
tine CALCPHOT to determine the corresponding broad-band 
magnitudes to explain the observed counts. The result, cor- 
rected for the small 12% renormalization found by Landsman 
(1997), was V = 28.80, Rc = 27.69 and Ic = 26.54. Although the 
uncertainty in the measured counts is only 10%, the consider- 
able dependence on the shape of the assumed spectrum, and the 
difficulty of absolutely calibrating the wide STIS passband (see 
below), leads us to adopt an uncertainty of 0.4 magnitudes for 
each band. 

To check the STIS result we tried an alternative calibra- 
tion using our WFPC2 photometry for the two brightest field 
stars, but still relying on SYNPHOT to calculate the appropriate 
colour correction, and found a STIS magnitude of V = 28.58. 
This discrepancy is larger than the formal errors, and is proba- 
bly indicative of the difficulty of precisely calibrating the very 
wide unfiltered STIS CCD passband. 

We collected and rereduced the near-infrared J- and K-band 
images of Soifer et al. (1997), obtained on March 30.2 and 31.2, 
1997 UT with NIRC on the Keck-1 10-m telescope, and the J-, 
H- and K -band images of Klose, Stecklum and Tuffs (1997), 
obtained with the Calar Alto 3.5-m telescope on March 17.8, 
1997 UT using the near-infrared camera MAGIC. The infrared 
frames were reduced by first removing bad pixels and combin- 
ing at least five frames around each object image to obtain a 
sky image. This sky image was then subtracted after scaling it 
to the object image level. 

A source was detected in both bands J and K in the Keck- 
1 observations. Relative to the nearby star (star H) we found 
J = 23.27±0. 15 (OT + GAL; sum of the March 30.3 and March 
31.2, 1997 UT data) and K = 22.85 ±0.25 (March 30.2, 1997 
UT) within a large 0."75 aperture radius. The J-band magnitude 
is consistent with the earlier determination (Soifer et al. 1997), 
but our K-band determination is 0.8 mag fainter Fruchter et al. 
(1999) detected the host galaxy at 7*: = 22.8 ± 0.3 and so our 
result is consistent with detection of the host galaxy only. 

We did not detect the source in the Calar Alto J, H and 
K images and determined the following limiting magnitudes: 
J > 21.45, H > 20.02 and K > 20.13 (3 a). We confirm the 
non detections in H and K by Klose et al. (1997), but do not 
confirm their detection of the source in J^"*. 

2.3. The host galaxy 

The host galaxy of GRB 970228 is clearly detected in the 
HST WFPC2 images, but at low signal-to-noise. Working with 
cosmic-ray and hot-pixel cleaned images, we first subtracted 

warm pixel at the position of the source that was not obvious in the bad pixel 



The J-band magnitude reported by Klose et al. (1997) is affected by a transient 
mask. 
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point-spread functions positioned and scaled to remove the im- 
age of the OT from each. For the two F606W images we then 
performed aperture photometry in apertures of increasing size, 
taking a sky estimate from the surrounding area which was free 
of other objects. We found that the counts reached an asymp- 
tote, within the noise, by a radius of 25 pixels We 
calibrated this using the procedure of Holtzman (1995b), find- 
ing an average Veoe = 25.5 ± 0.3, where the error is dominated 
by the photon and read noise. 

In the F814W image the host galaxy is even fainter, and we 
estimated its color by comparing the counts within a 7 pixel 
radius aperture with the counts in the same aperture on the 
F606W image. Thus we found a color Vfeoe-Zgu ~ 1.1 which is, 
correcting for the Galactic extinction, reasonable for a galaxy at 
this magnitude (Small et al. 1995). This could then be used to 
estimate color corrections to the standard systems, and we con- 
cluded V = 25.75 ± 0.35 and = 24.7 ± 0.5. 

The host galaxy is much better detected on the STIS image, 
although we must assume a spectral shape in order to estimate 
broad-band magnitudes. Using CALCPHOT we found a good 
match to the WFPC2 colours with a suitably redshifted and red- 
dened Sc galaxy spectrum'^. Based on this, we infer V = 25.77, 
Rc = 25.22 and = 24.73 from the STIS data, which is in good 
agreement with the other determinations. We adopt an error of 
0.2 mag to account for the systematic uncertainty introduced by 
the necessity of assuming the spectral shape. 

We did not detect the host galaxy of GRB 970228 in the ESQ 
VLT-UTl Science Verification images and determined the 3- 
cr limit for detection using an aperture radius of 1.5 times the 
seeing, which corresponds to ^ l."5. The results are shown in 
Table nl The V and Rc-band upper limits for the host galaxy 
show that the host galaxy is fainter than the earlier determina- 
tion (Sahu et al. 1997), but consistent with the result of Fruchter 
et al. (1999), Castander and Lamb (1998a) and with our esti- 
mate from the WFPC2 images. 

3. THE LIGHT CURVES 

In the interpolations to the Rc band (see Fig. |l]) we have 
assumed that the spectra of both the point source and the host 
galaxy, are smooth (i.e., not dominated by emission lines). We 
have used the relation between the color indices V -Rc and 
V -Ic given by The et al. (1984) for late-type stars; for bluer 
stars we have inferred this relation from the tables given by 
Johnson (1966) for main-sequence stars and the color transfor- 
mations to the Cousins VRI system given by Bessel (1976). 
We have tested the validity of these color-color relations from 
numerical integrations of power law flux distributions and of 
Planck functions, and conclude that if the flux distributions of 
the OT and the host galaxy are smooth, the uncertainty in the in- 
terpolated Rc magnitude is unlikely to exceed 0. 1 mag (Galama 
et al. 1997). For the HST September 4, 1997 observation we 
have assumed that the colors of the OT remained constant dur- 
ing the late-time decay (i.e., taking the observed V -1^ = 2.24 
from t he H ST March 26 and April 7, 1997 observations; see 
Sect. 2.2). Finally, we corrected for the contribution of the 
host galaxy emission (V = 25.77 ± 0.20, R^ = 25.22 ± 0.20 and 
/c = 24.73 ± 0.20) and obtained the V-, Rc-, and Ic-band Hght 
curves of the OT (Fig. |l]). 

We have fitted power laws, F oc f to the Ught curves (ignor- 



ing upper limits) and found: Fy = {11.7 ±1 .0)1^ 



= 0.2 with 2 d.o.f), Fr^ =12.7;|;2f-i i4±004. 



1.32±0.04 
,.2 _ 



(x-" = 9.7 with 8 



d.o.f), and =15.2!^;if7-^^*""VJy (x = 4.5 with 2 d.o.f), 
where fd is the time in days. These fits are indicated with dotted 
lines in Fig. |l[ We also fitted a power law to the Rc-band light 
curve at intermediate times (between 3 and 50 days) and found 

Fr^ =3.4;2 jf^'"^^° 'VJy ix^ = 1-3 with 4 d.o.f). The fit is in- 
dicated with a thick line in Fig. |l]. For the late-time detections 
(after 20 days) we found: Fy = (llJ± 1.3)fJ^"=^°-2'MJy (x^ 
= 0.2 with 1 d.o.f), Ft^ =21.5!23f-i 27±o.20^jy ^^2 ^ q 3 ^jjj^ j 

d.o.f) and Fi^ =50!^^f;['-^^=^'' 'Vjy (X^ = 0.9 with 1 d.o.f), where 
fd is the time in days. These fits are indicated with dashed lines 
in Fig. 

4. OPTICAL/NEAR-INFRARED TO X-RAY SPECTRA 

For three epochs we have reconstructed the spectral flux dis- 
tribution of the OT, corresponding to: (i) 1997 February 28.99 
UT (V, Rc and Ic data from Table and 2-10 keV flux from 
Costa et al. 1997), (ii) March 9.89 UT (Rc band measurement 
from Table |l| and ROSAT 0.1-2.4 keV unabsorbed flux from 
Frontera et al. 1998), and (iii) March 30.8 UT (See Fig. ||). 
We subtracted the host galaxy flux from ground-based measure- 
ments, corrected the OT fluxes for Galactic foreground absorp- 
tion. Ay = 0.78 ± 0.12 (from Schlegel, Finkbeiner and Davis 
1998; but see also the discussion by Gonzalez, Fruchter, and 
Dirsch 1999, and Castander and Lamb 1998b), and brought 
measurements to the same epoch using the power-law fits of 
Sect. ^ We note that the spectral flux distribution on March 
30.8 UT (Fig. |) is different from that of Reichart (1999) be- 
cause in our reanalysis the optical transient was not detected in 
the K band. The resulting spectral flux distribution therefore 
supersedes that of Reichart (1999) who obtained the K-band 
magnitudes of the OT and the host from the literature. 

We fitted the optical to X-ray spectral flux distribution of the 
first epoch (February 28.99, 1997 UT) with a power law and 
an exponential optical extinction law, /v oc v^e'^ , where we 
assumed that the extinction optical depth, r oc i/. The data did 
not allow determination of the host galaxy extinction; the fit 
provided a negative extinction {Ay < 0.50; 90 % confidence) 
and was subsequently fixed at zero. We found a spectral slope, 
/? = -0.780 ±0.022 (reduced x-squared, =4/2). For compar- 
ison, the X-ray power-law photon index observed ^12 hours 
after the burst (2-10 keV; Costa et al. 1997) corresponds to a 
spectral slope f3 = -1.1 ±0.3. 

5. DISCUSSION 

Most X-ray and optical/infrared afterglows display a power 
law temporal decay; this is a prediction for relativistic blast- 
wave models of GRB afterglows, in which a relativistically ex- 
panding shock front, caused by an energetic explosion in a cen- 
tral compact region, sweeps up the surrounding medium and ac- 
celerates electrons in a strong synchrotron emitting shock (e.g., 
Meszaros and Rees 1997). Assuming a single power law and 
ignoring upper limits for the moment, the V-, Rc- and Ic-band 
light curves (Fig. [l|) indicate that the power law decay rate is 
an increasing function of frequency; for the V, Rc and Ic band 
the power-law decay index a = -1.32 ±0.04, -1.14 ±0.04, and 
-0.92 ± 0.06, respectively. This implies that as a function of 
time the OT becomes redder and redder. Such a relation has not 
been observed previously, and is not what is expected in stan- 
dard relativistic blast-wave models, where the power law index 
is, to first order, independent of frequency. 



'^From the Kinney-Calzetti atlas, see the Synphot User Guide 1998, Space Telescope Science Institute. 
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Taking the upper limits into account, we find that the March 
4.87, 1997 UT NOT V-band non-detection is significantly be- 
low the fit. Also, the Rc-band light curve suggests that the op- 
tical emission is weaker around that time (between 3 and 10 
days) than expected from a single power law. The V and Rc 
data suggest that the initial decline was faster than the decline 
at intermediate times. From the V-band February 28.99 and 
March 4.87, 1997 UT data we find that the early-time decay 
was steeper than a = -1.4 {3a), while at intermediate times 
(around ^ 10 days) it was a = -0.73 ±0.17 (Rc). This struc- 
ture in the light curve was recognized early on by Galama et al. 
(1997), who showed that the early decay of the light curves (be- 
fore March 6, 1997) was faster than that at later times (between 
March 6 and April 7, 1997). The earlier analysis of the March 
26 and April 7, 1997 HST WFPC2 images yielded brighter es- 
timates of the host galaxy magnitudes (Sahu et al. 1997) than 
our reanalysis (see also Fruchter et al. 1999; Castander and 
Lamb 1998a). Because of this, the correction to the ground- 
based data for the host galaxy emission applied in this work is 
smaller than that of Galama et al. (1997), and thereby the effect 
is somewhat reduced, but it is still present (see also Reichart 
1999) and consistent with the earlier result. 

The late-time light curves show that the rate of decay in- 
creased again after ^ 35 days. The late-time decays in V, R^ and 
Ic are consistent with being frequency independent (i.e., simi- 
lar decays in V, Rc and Ic), with a power-law index a = -1.2. 
The back extrapolation of the late-time fit in the V-band pre- 
dicts quite well the earliest V-band data, but this is not the case 
for the Rc- and Ic-band light curves. We conclude that a single 
power-law is not a good representation of the light curves; the 
overall shape of the lightcurves is that of an initial fast decay 
(before ^ March 6, 1997) followed by a 'plateaux' (slow de- 
cay; between ^ March 6 and ^ April 7, 1997) and finally the 
light curves resume a fast decay (after ^ April 7, 1997). 

For synchrotron radiation, by a power-law energy distribu- 
tion of electrons, at low frequencies, and for an adiabatic evo- 
lution of the blast wave we can distinguish the following two 
cases, each of which has its own relation between the spec- 
tral slope /3 and the power law decay index, a (Sari, Piran 
and Narayan 1998): (i) both the peak frequency and the 
cooling frequency are below the optical waveband. Then 
a = 3/3/2+1/2 = -0.67 ± 0.03 (taking the measured spectral 
slope, (3 = -0.78 on February 28.99, 1997 UT) (ii) has 
passed the optical waveband, but has not yet. In that case 
a = 3/3/2 = -1 . 17± 0.03. The temporal decay values measured 
at early times are in reasonable agreement with case (ii), but 
certainly not with case (i), suggesting that the peak frequency 

is below the optical passband and that the cooling frequency 
t^c is near or above the X-ray (2-10 keV) passband; such a high 
cooling frequency appears to be the fairly general case in GRB 
afterglows (e.g.. Bloom et al. 1998; Vreeswijk et al. 1999; 
Galama et al. 1999). The fit to the optical to X-ray spec- 
trum of February 28.99, 1997 UT suggests that extinction at 
the host is not very significant (Ay < 0.5; corresponding to a 
restframe extinction Ay < 0.3). At the second epoch we mea- 
sured a spectral slope, (3 = -0.52 ± 0.08, on March 9.89, 1997 
UT. For case (ii) this would correspond to a power-law decay of 
a = -0.78 ±0.12, consistent with a shallower decay at interme- 
diate times. We thus find that the observations of the early-time 



afterglow are roughly consistent with the predictions of sim- 
ple relativistic blast-wave models (see also Wijers, Rees and 
Meszaros 1997; Waxman 1997; Reichart 1997), except that we 
find evidence for a flattening of the decay rate around ^ 10 
days. 

As mentioned above, the frequency dependence of the rate of 
decay of the light curves corresponds to a reddening of the opti- 
cal afterglow with time (V-Zc = 0.50 ±0.23 at February 28.99, 
1997 UT and V-k = 2.24 ±0.14 at March 26 and April 7, 

1997) . We are confident about this large change in color as the 
HST calibration of three field stars agrees well with that of the 
ground-based calibration. Assuming a power-law spectrum, the 
optical spectral slopes (Galactic extinction corrected) from the 
intermediate-time V and Ic HST observations are (3 = -3.8 ±0.7 
(March 26, 1997) and /3 = -4.0 ±0.8 (April 7, 1997). These 
steep spectra cannot be explained by extinction at the host 
galaxy; we find a staggering Ay = 4.5 ± 1.0 (fitting as before for 
an exponential optical extinction law, /v « t^'^e"^, with t (x i/, 
and fixing the spectral slope at /? = -0.78), which would imply 
that the extinction increased with time from negligible (Ay < 
0.5) to very significant (Ay = 4.5 ± 1.0); this seems quite un- 
likely. A steepening of the spectral slope may be expected from 
a number of spectral transitions: e.g, the passage of the cooling 
frequency or the transition to a Sedov-Taylor non-relativistic 
expansion (e.g.. Sari et al. 1998; Waxman, Kulkarni and Frail 

1998) . However, such spectral transitions will be accompanied 
with an increase in the rate of decay of the light, which we do 
not observe here. The discrepancy is even larger if one consid- 
ers the spectral flux distribution around the third epoch (March 
30.8, 1997 UT; Fig. Eh; the spectrum shows an unexpected 
turnover around 3 x 10 Hz (Fruchter et al. 1999). This break 
plus the very red spectrum in V and Ic in the intermediate-time 
optical flux distribution is not consistent with the predictions of 
simple relativistic blast-wave models (see also Reichart 1999; 
Fruchter et al. 1999). 

What then could explain these observations? Bloom et al. 
(1999) observed a similarly red spectrum for the late-time emis- 
sion from GRB 980326. The authors argue that this is due to a 
supernova (SN) that dominated the light at late times. Reichart 
(1999) suggests that this is also the case for GRB 970228. In 
this hypothesis the early light curve is dominated by a (power- 
law decaying) GRB afterglow (produced by the relativistic blast 
wave), while at late times the light curve contains a significant 
contribution from an underlying supernova. Such would be 
the natural outcome of the 'collapsar/microquasar' models in 
which GRBs arise in jets that are formed in the core collapse to 
a black hole of massive stars (Woosley 1993; Paczyhski 1998). 
Such massive stars will, prior to collapse, have lost their hydro- 
gen envelopes and so we expect the supernova to be of type I^/c 
(Woosley 1993; Woosley, Eastman and Schmidt 1999; Bloom 
et al. 1999). 

Following Bloom et al. (1999) and Reichart (1999) we con- 
structed simulated supernova light curves and spectral flux dis- 
tributions as follows. First we combined three sets of Galactic 
extinction corrected U-, B-, V-, Rc- and Ic-band observations of 
the well-observed, unusual type-Ic SN 1998bw (Galama et al. 
1998; McKenzie and Schaefer 1999; Vreeswijk et al. 1999), 
which is likely to be associated with GRB 980425 (Galama et 
al. 1998; Kulkarni et al. 1998)"'. We then redshifted the 



'^The association between SN 1998bw and GRB 980425 lias been considered uncertain because the behavior of a another source in the eiTor box of GRB 980425, 
the fading X-ray source ISAX J1935. 3-5252, resembled that of a typical GRB afterglow. However, reanalysis of NFI BeppoSAX observations of the GRB 980425 
field by Pian et al. (1999) shows that its behavior is unlike that of previously observed GRB afterglows. SN 1998bw remains not typical for GRB afterglows because 
of its low redshift, z = 0.0085. 
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lightcurves to the redshift of GRB 970228, z = 0.695 (Djorgov- 
ski et al. 1999), which consists of wavelength shifting, time- 
profile stretching, and dimming (we assumed a cosmology with 
H(, = 70 km s"' Mpc"' and JIq = 0.3). To obtain the correspond- 
ing redshifted light curves in the V, Rc and Ic bands we assumed 
a power-law spectral shape to interpolate (and extrapolate) the 
redshifted SN 1998bw light curves to the central wavelengths of 
the V, Rc and Ic bands. (Note that the redshifted central wave- 
length of the U band is near V and Rc, and for the B band is near 
Ic.) Finally we corrected the SN light curves for the Galactic 
extinction in the direction of GRB 970228 {Ay = 0.78 ± 0.12). 
These light curves are shown in Fig. ||. 

The supernova spectral flux distribution corresponding to the 
epoch of March 30.8, 1997 UT (30.7 days after GRB 970228) 
is shown in Fig. ^ (here we did the opposite: we corrected 
the GRB 970228 spectral flux distribution for Galactic extinc- 
tion). The resemblance of the GRB 970228 and the redshifted 
SN 1998bw spectral flux distributions is remarkable (see also 
Reichart 1999). 

At early times the Rc-band light curve (see Fig. |[) will not be 
very much dominated by SN emission. If we attribute the early 
part (f < 8 days) of the Rc-band light curve to GRB afterglow 
and fit a power-law temporal decay to the data we find a tempo- 
ral decay index a = -1.46 ±0.16 (x^ =1.5 with 3 d.o.f). This 
value of the decay constant is consistent with that in X rays, 
a = -1.33![S [^ during the first 4 days (2-10 keV; Costa et al. 
1997) and a = -1.50!j;:| up to - 10 days (0. 1-2.4 keV; Frontera 
et al. 1998). It is also consistent with the earlier determination 
(Galama et al. 1997; see also Reichart 1999). 

Next we took a single value for the power-law decay in V, 
Rc and Ic, computed the sums of these power-law decays plus 
supernova light curves, and fitted for the power-law decay in- 
dex, a, by minimizing the x-squared of the V, Rc, Ic light 
curves and the resulting fit. But first, the assumption of a con- 
stant color that was used to derive a V-band magnitude from 
the September 1997 unfiltered STIS observation is no longer 
valid as in this hypothesis at late times the SN emission will 
dominate the light and so we expect the light to have reddened 
(see Fig. Thus, in order to convert the observed counts 
into broadband magnitudes we assumed that the spectral en- 
ergy distribution was similar to that of SN 1998bw, redshifted 
to z = 0.695. Our procedure was as follows: (i) we used the red- 
shifted U, B, V, Rc, Ic spectral-flux distribution of SN 1998bw 
(on day 111, corresponding to the STIS epoch at 188 days), (ii) 
we corrected for Galactic extinction. Ay = 0.78 ±0.12, interpo- 
lating to the redshifted central wavelengths of the U, B, V, Rc, 
Ic bands using the Galactic extinction curve of Cardelli, Clay- 
ton and Mathis (1989) to obtain the spectral flux distribution as 
it would be observed in the direction of GRB 970228, and (iii) 



we used this spectral flux distribution in the SYNPHOT rou- 
tine CALCPHOT, and applied the Landsman (1997) correction, 
to determine corresponding broad-band magnitudes, V = 28.56, 
7?c = 27.42 and 4 = 26 A3. Again, we adopt an uncertainty of 
0.4 mag to make some allowance for the sensitivity of the re- 
sult to the assumed spectral shape. Finally we minimized the 
X-squared and found a = -1.73![!:°^ (X^ = 17.5 with 14 d.o.f; 
we note that the x-squared is not strictly correct since a few 
of the data points are not independent). The fit is excellent in 
view of the uncertainties in the model. Similar results were 
presented by Reichart (1999). The fit is somewhat better than 
that presented by Reichart (1999), presumably because of our 
consistent reanalysis of all the original images, and, in partic- 
ular, because our STIS magnitudes were calculated using the 
knowledge of the spectral shape of the template supernova. 

The model has very few assumptions: the power law index a 
is the only free parameter and we did not allow the peak lumi- 
nosity of the SN to be a free parameter but used the SN 1998bw 
observations as they are. SN 1998bw is, of course, an unusual 
example of the class, being very bright and rapidly expanding, 
but its coincidence with GRB 980425 makes it the only tem- 
plate we can reasonably choose for now. Also, from the 'col- 
lapsar' model we expect a GRB to be accompanied by a su- 
pernova of type Ib/c. In any case, the light curves (Fig. |3]) as 
well as the intermediate-time spectral flux distribution (Fig. ^ 
and intermediate-time colors are remarkably well explained by 
the sum of a power-law plus supernova light curve. It can ex- 
plain the apparent transition from a rapid to a slow and again a 
rapid decline of the light curves, is in good agreement with the 
upper limits in the light curves and it can account for the un- 
usual intermediate-time red color V -1^ = 2.24 and the peaked 
spectral flux distribution on March 30.8, 1997 UT (see also Re- 
ichart 1999). These features cannot be accounted for by simple 
relativistic blast-wave models. We thus find the data to be in 
agreement with the hypothesis of Reichart (1999) that a super- 
nova dominated the light at late times. Together with the evi- 
dence for GRB 980326 (Bloom et al. 1999) and GRB 980425 
(Galama et al. 1998) this is further support for the idea that 
at least some GRBs are associated with a possibly rare type of 
supernova. 
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NASA grant NAG 5-2560. 



REFERENCES 



Baggett, S., Casertano, S., Gonzaga, S., Ritchie, C. 1997, Instrument Science 

Report WFPC2 97-10 
Bessel, M.S. 1976, PASP, 88, 557 
Bloom, J.S. et al. 1998, ApJ, 508, L21 

Cardelli, J.A., Clayton, G.C. and Mathis, J.S. 1989, ApJ, 345, 245 
Casali, M.M. and Hawarden, T.G. 1992, JCMT-UKIRT Newsletter, No. 3, 33 
Casertano, S. and Mutchler, M. 1998, Instrument Scievire Rppnrt WFPr]^ 98-02 
Castander, F.J. and Lamb, D.Q. 1998a, in press, ApJ, |pstrn-ph/9807l9'^ 
Castander, F.J. and Lamb, D.Q. 1998b, in press, ApJ,]astro-phyy8U / 19* 

Costa, E. et al. 1997, Nature, 387, 783 

Djorgovski, S.G., et al. 1997, lAU Circ, 6732 
Djorgovski, S.G. et al. 1999, GCN Circular, No. 289 
Frontera, F. et al. 1998, A&A, 334, 69 
Fruchter, A.S. et al. 1997, lAU Circ, 6747 



Fruchter, A.S. et al. 1999, ApJ, 516, 683 
Galama, T.J. et al. 1997, Nature, 387, 479 
Galama, T.J. et al. 1998, Nature, 395, 670 
Galama, T.J. et al. 1999, Nature, 398, 394 

Gonzalez, R.A., Fruchter, A.S., and Dirsch, B. 1999, ApJ, 515, 69 
Groot, RJ. et al. 1997a, lAU Circ, 6584 
Groot, RJ. et al. 1997b, lAU Circ, 6588 
Guarnieri, A. et al. 1997, A&A, 328, L13 
Holtzman, J. et al. 1995a, PASP, 107, 156 

Holtzman, J. A., Burtows, C.J., Casertano, S., Hester, J.J., Trauger, J.T., Watson, 

A.M., Worthey, G., 1995b, PASP, 107, 1065 
Johnson, H.L. 1966, ARA&A, 4, 191 
Klose, S., Stecklum, B., Tuffs, R. 1997, lAU Circ, 661 1 
Kulkami, S.R. et al. 1998, Nature, 395, 663 



7 



Landolt, A.U. 1992, AJ, 104, 340 

Landsman, W. 1997, STIS CCD image mode sensitivity. Internal STIS 

Calibration Document. 
Margon, B., Deutsch, E.W., Lamb, D.Q., Castander, F.J. 1997, lAU Circ, 6618 
McKenzie, E.H. and Schaefer, B.E. 1999 
Meszaros, P. and Rees, M.J. 1997, ApJ, 476, 232 
Metzger, M.R. et al. 1997a, lAU Circ, 6588 
Metzger, M.R. et al. 1997b, lAU Circ, 6631 
Paczynski, B. 1998, ApJ, 494, L45 
Pedichini et al. 1997, A&A, 327. L32 
Plan, E. et al. 1999, A&AS, in press 

Reichart, D.E. 1997, ApJ, 485, L57, . 

Reichart, D.E. 1999, ApJ, in press, astro-ph/990607g 
Sahu, K. et al. 1997, Nature, 387, 47B 



Sari, R., Piran, T., and Narayan, R. 1998, ApJ, 497, L17 
Schlegel, D.J., Finkbeiner. D.P and Davis, M. 1998, ApJ500, 525 
Soifer, B. et al. 1997, lAU Circ, 6619 

Small, I., Hogg, D.W., Yan, L. and Cohen, J.G. 1995, ApJ, 449, L105 
The, PS., Steenman, H., Alcaino, G., 1984, A&A, 132, 385 
Van Paradijs, J. et al. 1997, Nature, 386, 686 
Vreeswijk, P.M. et al. 1999, in press, ApJ 
Waxman, E. 1997, ApJ, 485, L5 

Waxman, E., Kulkarni. S.R. and Frail. D.A. 1998, ApJ, 497, 288 
Whitmore, B. and Heyer, I. 1997, Instrument Science Report WFPC2 97-08 
Whitmore, B. 1998, Technical Instrument Report WFPC2 98-01 
Wijers, R.A.M.J., Rees, M.J., and Meszaros. P 1997, MNRAS, 288, L51 
Woosley, S.E. 1993, ApJ, 405, 273 

Woosley, S.E., Eastman, R.G. and Schmidt, B.P 1999, ApJ, 516, 788 



Table 1 

Summary of optical observations. Upper limits are three sigma. Abbreviations: OT, optical transient; GAL, host galaxy; RAO, 
Rome Astrophysical Observatory; BUT, Bologna University Telescope; WHT, William Hershell Telescope; APO, Apache Point 
Observatory; NOT, Nordic Optical Telescope; INT, Isaac Newton Telescope; CA, Calar Alto 3.5-m telescope; HST, Hubble Space 
Telescope; P5m, Palomar 200-inch; VLT-UTl, Very Large Telescope Unit I. 



Date (UT) 


Telescope 


Exp. time (sec) 


Magnitude 


Remarks 


Feb. 28.81 1997 


RAO 




^c = 20.0 ±0.5 


Unfiltered OTh-GAL 


Feb. 28.83 1997 


BUT 




/?c = 20.58 ± 0.28 


OTh-GAL 


Feb. 28.99 1997 


WHT 




y = 21.01±0.10 


OTh-GAL 


Feb. 28.99 1997 


WHT 




/c = 20.51 ±0.21 


OTh-GAL 


Mar 03.10 1997 


APO 




Bj = 24.5 ±0.7 


OT+GAL 


Mar 03.90 1997 


NOT 




> 22.70 ±0.20 


Unfiltered OT+GAL 


Mar 04.86 1997 


NOT 




/?e = 23.28 ±0.40 


Unfiltered OT+GAL 


Mar 04.87 1997 


NOT 




V > 23.96±0.12 


OT+GAL 


Mar 04.89 1997 


NOT 




fi> 23.87±0.10 


OT+GAL 


Mar 06.32 1997 


Keck 




/?c = 23.60 ±0.22 


OT+GAL 


Mar 08.86 1997 


INT 




y> 24.00 ±0.10 


OT+GAL 


Mar 08.88 1997 


WHT 




/c > 22.68 ±0.10 


OT+GAL 


Mar 08.89 1997 


INT 




i^c > 23.35 ±0.08 


OT+GAL 


Mar 09.85 1997 


INT 




B = 25.43 ±0.31 


OT+GAL 


Mar 09.89 1997 


INT 




i?c = 23.75 ±0.28 


OT+GAL 


Mar 13.00 1997 


NTT 




^c = 24.36 ±0.15 


OT+GAL 


Mar 17.84 1997 


CA 




7>21.45±0.15 


OT+GAL 


Mar 17.83 1997 


CA 




//> 20.02 ±0.15 


OT+GAL 


Mar 17.81 1997 


CA 




k' > 20.13 ±0.15 


OT+GAL 


Mar 26.38 1997 


HST 




y = 25.98 ±0.22 


OT 


Mar 26.47 1997 


HST 




/c= 23.76 ±0.11 


OT 


Mar. 30.27 1997 


Keck 




= 22.85 ±0.25 


OT+GAL 


Mar. 30.30 H- 31.26 1997 


Keck 




7 = 23.27±0.15 


OT+GAL 


Apr. 05.76 1997 


Keck 




= 24.97± 0.25 


OT+GAL 


Apr. 07.22 1997 


HST 




y = 26.36±0.26 


OT 


Apr. 07.30 1997 


HST 




4 = 24.10 ±0.15 


OT 


Mar. 26 + Apr. 7 1997 


HST 




y = 25.75 ±0.35 


GAL 


Mar. 26 + Apr. 7 1997 


HST 




/c = 24.7 ±0.5 


GAL 


Sep. 04 1997 


P5m 




/? = 25.5 ±0.5 


GAL 


Sep. 04.71 1997 


HST 




y = 28.80±0.40 


OT 


Sep. 04.71 1997 


HST 




i?c = 27.69 ±0.40 


OT 


Sep. 04.71 1997 


HST 




/c= 26.54 ±0.40 


OT 


Sep. 04.71 1997 


HST 




y = 25.77 ±0.20 


GAL 


Sep. 04.71 1997 


HST 




7^0 = 25.22 ±0.20 


GAL 


Sep. 04.71 1997 


HST 




7^=24.73 ±0.20 


GAL 


Sep. 01.35 1998 


VLT-UTl 


1800 


B> 26.08 ±0.08 


GAL 


Sep. 01.35 1998 


VLT-UTl 


1200 


y> 25.33 ±0.05 


GAL 


Sep. 01.35 1998 


VLT-UTl 


1500 


/?c > 25.29 ±0.06 


GAL 
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Table 2 

B, V, Rc, Ic magnitudes, associated random errors AB, AV, ARc, Ale, and the coordinates of the reference stars. We estimate the 
absolute calibration uncertainty to be 0.05 magnitudes in B, V, and Rc and 0.10 magnitudes in Ig. 



Star 




A 


B 


C 


D 


E 


F 


G 


H 


B 




21.693 


19.614 


23.161 


23.654 


22.296 


23.976 






AB 




0.022 


0.010 


0.042 


0.060 


0.016 


0.073 






V 




20.345 


18.793 


21.782 


21.972 


20.648 


22.137 


23.237 


23.003 


AV 




0.011 


0.010 


0.017 


0.021 


0.010 


0.023 


0.064 


0.044 


Rc 




19.646 


18.413 


20.874 


21.119 


19.640 


20.889 


23.008 


21.906 


ARe 




0.040 


0.040 


0.025 


0.032 


0.010 


0.035 


0.103 


0.067 


Ic 




18.787 


17.758 


20.099 


20.501 


18.665 


19.279 




20.923 


Ale 




0.034 


0.034 


0.045 


0.062 


0.034 


0.040 




0.080 


RA (J2000) 


5:01: 


48.9 


49.3 


47.8 


48.7 


45.9 


45.6 


47.4 


46.4 


Dec (J2000) 


11: 


45:59 


45:43 


46:56 


47:01 


47:12 


47:24 


46:56 


46:54 
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Fig. 1. — The V-, Rc- and Ig-band lightcurves of GRB 970228 (magnitude versus time). Rc-band data obtained from interpolations 
between V and Ic are indicated with open symbols (see the discussion for details). Indicated are power-law fits, /v oc f", to the 
detections (ignoring upper Umits). Fits to all points are indicated by the dotted lines and fits to the late- time lightcurve (after 20 days) 
are indicated by the dashed lines (for the V band they overlap). Also indicated is a fit to the Rc-band Ught curve at intermediate times 
(between 3 and 50 days; thick Une). 
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Spectral flux distribution March 30.8, 1997 UT 
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Fig. 2. — The optical/near-infrared spectral flux distribution of GRB 970228 at March 30.8, 1997 UT (• and upper-limit arrow); 
from left to right: K-band detection (Table [l|) translated into an upper limit because it is consistent with detection of the host 
galaxy only (detection magnitude minus 3 a), J-band detection (Tablell|) translated into an upper limit due to unknown J-band host 
magnitude (detection magnitude minus 3 <t), V and Ic HST measurements from Table [l] Also shown is the spectral flux distribution 
of SN 1998bw redshifted to the redshift of GRB 970228 (o). Note that the spectral flux distribution is different from that of Reichart 
(1999) (See the discussion in Sect. ^. 



lifiht curves of GRB 970228 




1 10 100 1 10 100 1 10 100 



Days after Feb 28.1236 UT 

Fig. 3. — The V-, R^-, and Ic-band lightcurves of GRB 970228 (fluxes versus time). The dotted curves indicate power-law decays 
with a = -1 .73, and redshifted SN 1998bw light curves. The thick line is the resulting sum of SN and power-law decay light curves. 



